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TIME  RESOLVED  MEASUREMENTS  OF  LARGE  AMPLITUDE  VELOCITY 
FLUCTUATIONS  WITH  A  THREE-SENSOR  HOT-HIRE  PROBE 


W.  W.  Willmarth,  T.  Wei  and  K.  Madnia 
Department  of  Aerospace  Engineering 
j  University  of  Michigan,  Ann  Arbor,  Michigan  48109 

Calibration  and  data  reduction  procedures  are  reported  for  a 
three  sensor  hot-wire  probe  (TSI  model  1298)  used  to  determine 
all  three  velocity  components  in  large  amplitude  turbulence  at  a 
speed  of  65  ft./ sec.  Hotrwire  sensors  were  necessary  because 
hot-film  sensors  did  not  produce  a  steady,  predictable  signal 
for  a  large  variation  of  the  angle  between  the  sensor  axis  and 
the  stream.  A  digital  computer  was  used  to  obtain  an  iterative 
solution  of  the  equations  relating  the  sensor  signals  to  the 
velocity  components.  Reasonably  accurate  measurements  of  the 
time  resolved  velocity  components  were  made  for  transverse 
velocities  less  than  approximately  2/10  of  the  stream  velocity. 
For  larger  transverse  velocities,  the  instability  of  the 
iterative  solution  in  conjunction  with  the  occurrence  of  multiple 
solutions  presents  a  serious  problem.  Improved  methods  of 
iteration  and  additional  information  beyond  that  provided  by  the 
three  sensor  signals  appears  to  be  necessary  to  determine  the 
correct  solution  among  a  maximum  of  four  possible  solutions  with 
flow  downstream  along  the  probe  axis. 


INTRODUCTION 


to: 


This  paper  is  a  study  of  the  problem  of  making  time  resolved 
measurements  of  all  three  velocity  components  of  a  highly 
turbulent  flow  using  a  three-sensor  hot-wire  probe.  The 
literature  describing  the  use  of  hot-wires  and  other  types  of 
heated  sensors  for  turbulence  measurements  is  extensive.  The 
recent  review  articles  by  Comte-Bel lot1  and  Blackwelder2  and  the 
book  by  Perry-*  provide  detailed  information  on  the  hot-wire 
method  of  measurement  and  an  introduction  to  the  literature. 
Unfortunately  the  problem  of  making  accurate,  time  resolved 
measurements  of  all  three  large  amplitude  velocity  fluctuations 
using  three  hot-wires  is  not  treated  in  detail  in  the  above 
review  articles. 

When  three  dimensional  velocity  fluctuations  are 
encountered,  the  response  of  a  single  sensor  is  a  function  of 
all  three  velocity  components.  One  must  make  simultaneous 
measurements  of  the  output  signals  from  at  least  three  sensors, 
in  order  to  obtain  the  data  required  to  determine  all  three 
turbulent  velocity  components.  If  the  turbulence  level  is  low  a 
linearized  approximate  solution  of  the  three-sensor  hot-wire 
response  equations  is  possible.  However,  a  linearized  solution 
is  not  adequate  for  large  amplitude  turbulent  velocity 
fluctuations.  Digital  computers  fitted  with  suitable  A/D 
converters  are  a  natural  choice  for  the  simultaneous  measurements 
and  computations  required  for  measurements  of  this  type. 

In  1974,  Lakshminarayana  and  Poncet4  reported  the  use  of  an 
X-wire  probe  and  a  single-wire  probe  to  measure  three  velocity 
components  in  a  turbulent  wake  flow  behind  a  turbomachinery 


rotor.  The  data  from  each  sensor  was  digitized  and  a  digital 
computer  was  used  to  linearize  the  sensor  signals  thus 
determining  the  the  velocity  normal  to  each  sensor.  The  computer 
was  then  used  to  solve  three  simultaneous  equations  relating  the 
velocity  normal  to  the  sensors  to  the  actual  velocity  components. 
Average  values  of  the  mean  and  fluctuating  velocity  components 
were  reported  and  estimates  of  the  accuracy  of  the  results  based 
on  statistical  measures  for  a  Gaussian  process  were  made.  The 
method  and  accuracy  of  the  solution  of  the  three  simultaneous 
equations  and  the  possible  occurrence  of  multiple  roots  were  not 
documented . 

In  1978,  Moffat,  Yavuzkurt,  and  Crawford^  developed  and 
used  a  custom  built  analog  computer  to  obtain  the  instantaneous 
velocity  components  in  real  time  using  a  DISA  triaxia 1-wire 
probe.  In  the  same  year,  Fabris^  reported  the  development  of  a 
special  four-wire  probe  (one  wire  measured  temperature)  and  the 
use  of  a  digital  computer  employing  the  Newton-Raphson  technique 
to  iteratively  solve  three  simultaneous  nonlinear  hot-wire 
response  equations  in  a  heated  turbulent  wake  flow. 

Recently,  1983,  Paulsen7  described  a  method  for  obtaining 
a  solution  for  the  three  velocity  components  with  a  digital 
computer  using  a  multiple,  modified  Newton  search  algorithm. 
Paulsen7  also  reported  the  use  of  this  method  of  solution  for 
turbulence  data  obtained  with  a  subminiature  triple  sensor  probe 
of  his  own  design.  The  accuracy  of  the  data  was  not  satisfactory 
because  prong  interference  was  a  serious  problem  (the  prong 
spacing  was  three  times  the  prong  tip  diameter)  but  Paulsen 


stated  that  improvements  of  the  probe  or  corrections  for  the 
prong  interference  will  be  made  in  the  future. 

Chang  and  Adrian®,  in  1984,  described  the  use  of  the  half 
interval  method  to  obtain  a  digital,  iterative  solution  for  the 
response  of  a  three  wire  probe  to  the  velocity  components  in  a 
high  intensity  turbulent  flow  around  the  potential  core  of  a 
turbulent  jet.  Jorgensen's®  cooling  law  was  used  for  the  wire 
response  and  the  probe  was  a  standard  X-wire  configuration  with 
an  adjacent  slanted  wire.  Multiple  solutions  of  the  response 
equation  were  encountered  and,  as  the  turbulence  level 
increased,  the  occurence  of  non-realizable  and  non-sol vable  data 
was  documented.  The  time  resolved  solutions  for  the  velocity 
components  were  not  presented  nor  was  it  possible  to  compare  the 
computed  velocity  components  with  the  instantaneous  velocity 
components  to  which  the  probe  was  exposed  in  the  turbulent  flow. 

Wallace10  has  reported  the  use  of  an  array  of  nine  hot-wire 
sensors  for  measurements  of  the  velocity  and  vorticity  components 
at  a  point  in  a  turbulent  flow.  As  might  be  expected,  the 
calibration  and  data  processing  procedures  were  complex. 
According  to  Wallace,  in  a  private  communication,  problems 
similar  to  those  encountered  by  Chang  and  Adrian  occurred  and  the 
measurements  had  to  be  restricted  to  low  amplitude  turbulent 
fluctuations  to  avoid  serious  errors  in  interpretation  of  the 
probe  signals. 


In  the  present  paper  the  development  of  a  three  sensor  probe 
capable  of  measurements  of  large  scale,  large  amplitude 
turbulent  velocity  fluctuations  is  described.  The  accuracy  and 
limitations  of  the  calibration  and  data  reduction  procedures  are 
studied.  The  results  of  time  resolved  measurements  with  the 
probe  are  directly  compared  with  the  actual  velocity  to  which  the 
probe  is  exposed.  The  probe  calibration  and  data  reduction 
methods  have  been  applied  to  the  measurement  of  turbulent 
velocity  fluctuations  of  large  scale  in  the  atmosphere,  see 
Willmarth1* . 

I.  EXPERIMENTAL  EQUIPMENT 

A.  Apparatus 

A  three  sensor  probe  (model  1298  manufactured  by  TSI)  was 
selected  for  measurements  of  large  scale,  large  amplitude 
velocity  fluctuations  in  the  turbulent  flow  encountered  by  a 
vehicle  driven  on  an  expressway.  The  probe  was  calibrated  by 
pitching  and  yawing  it  in  the  flow  in  the  test  section  of  the  61 
X  61  cm.  wind  tunnel  at  the  Aerospace  Engineering  Department  of 
the  University  of  Michigan.  A  description  of  this  wind  tunnel 
has  been  given  by  Uberoi1^.  Fig.  1  is  a  drawing  of  the  probe 
mounted  at  the  center  of  the  tunnel  on  two  model  aircraft  servo 
units  (Airtronics  model  94554)  which  are  capable  of  a  maximum 
angular  rotation  of  plus  or  minus  60  degrees.  The  pitch  and  yaw 
servo  signals  were  generated  by  manually  rotating  the  shafts  of 
two  potentiometers  of  a  servo  signal  generating  device.  The 
resulting  angular  orientation  of  the  servo  shafts,  ©p  ,  ©y  , 
(see  Fig.  1)  was  determined  from  the  signals  produced  by 


two  additional  potentiometers  (not  shown)  mounted  on  each  servo 
shaft. 


The  three  velocity  components  relative  to  the  probe,  as 
defined  in  Fig.  1,  are  given  by: 

u  =  U  ^  cos  cos  9y  (1) 

v  =  U ^  sin  0p  (2) 

w  «  cos^0psin©y  (3) 

where,  ©  ,  is  the  pitch  angle,  ©,  is  the  yaw  angle  and,  U  , 
P  7  w 

is  the  free  stream  velocity  along  the  axis  of  the  wind  tunnel. 
The  free  stream  velocity  was  measured  with  a  manometer  connected 

to  static  pressure  taps  in  the  settling  chamber  and  the  test 

section  of  the  wind  tunnel. 

The  probe  sensors  were  operated  in  the  constant  temperature 
mode  using  TS1  model  1054B  constant  temperature  linearized 
anemometers.  The  linearized  sensor  signals  along  with  the  servo 
shaft  pitch  and  yaw  potentiometer  signals  were  recorded  on  analog 
magnetic  tape  using  a  Honeywell  model  5600C  tape  recorder  which 
had  seven  frequency  modulated  record  and  reproduce  channels  with 
frequency  response*  from  zero  to  5  KHz.  The  recorded  signals, 
for  various  flow  speeds  and  probe  orientations,  were  reproduced 
and  digitized  using  a  12  bit  analog  to  digital  converter 
(Analogic  model  5800)  controlled  by  a  digital  computer  system 
(Data  General  NOVA  840).  The  computer  system  was  also  used  to 
perform  the  computations  described  later  in  the  paper. 


B.  Three  Sensor  Probe 

Fig.  2  is  a  sketch  of  the  three  sensor  probe  (TSI  model 
1298)  showing  the  velocity  components  u,  v,  and  w  relative  to 
the  probe.  Initially,  a  probe  fitted  with  hot-film  sensors 
(TSI-60)  was  purchased  since  sturdy  sensors  were  desired  for  the 
measurements  to  be  made  on  the  highway.  As  described  below,  the 
hot-film  sensors  were  inadequate  for  accurate  measurements  of 
large  amplitude  velocity  flucuations.  Hot-wire  sensors  (TSI  T- 

I. 5)  which  proved  to  be  adequate  for  the  measurements,  are 
shown  mounted  on  the  prongs  of  the  probe  and  are  numbered  from  1 
to  3.  The  probe  was  slightly  modified,  as  discussed  in  section 
(III  A),  to  provide  a  larger  spacing  between  the  sensors. 

II.  RESULTS  FOR  CALIBRATION  OF  PROBE  FITTED  WITH  HOT-FILM  SENSORS 

A  standard  TSI  model  1298  probe  fitted  with  hot-film  sensors 
(TSI-60)  was  mounted  on  the  servo  positioner  and  exposed  to  a 
uniform  flow  of  19.4  m/sec,  comparable  to  the  maximum  speed  to  be 
used  for  measurements  on  an  expressway.  Linearized  signals  from 
the  hot-film  sensors  and  the  signals  from  the  servo  pitch  and  yaw 
potentiometers  were  recorded  when  the  probe  was  pitched  with  zero 
yaw  and  yawed  with  zero  pitch.  It  was  observed  that  the  traces 
of  the  sensor  outputs,  plotted  on  an  x-y  plotter  as  functions  of 
pitch  and  yaw,  see  Fig.  3  and  4,  were  irregular.  The  trace  from 
sensor  (2)  when  the  probe  was  pitched  at  zero  yaw  and  from  sensor 
(3)  when  the  probe  was  yawed  at  zero  pitch  exhibit  the  greatest 
irregularity.  There  is  a  sudden  increase  in  the  linearized 
output  voltage  as  the  pitch  or  yaw  angle  increases. 

At  first  it  was  thought  that  these  irregularities  were 


caused  by  non  uniform  flow  in  the  wind  tunnel.  However,  further 
examination  of  the  sensor  signals  at  discrete  angles  of  pitch  and 
yaw  revealed  that  they  were  not  constant  and  contained  sinusoidal 
oscillations  at  various  frequencies  from  2  to  7  KHz.  These 
oscillations  appeared,  changed  and  disappeared  at  various  angles 
of  pitch  and  yaw.  The  Reynolds  number  was  approximately  218 
based  on  the  sensor  diameter  and  the  average  component  of  the 
19.4  m/sec.  velocity  normal  to  the  slanted  sensors.  It  is 
apparent  that  the  rate  of  heat  transfer  from  the  sensor  is  not  a 
predictable  function  of  the  angle  between  the  axis  of  the  sensor 
and  the  free  stream  velocity.  It  was  concluded  that  the  152yt(m 
diameter  hot-film  sensors  (TSI-60)  were  not  suitable  for  use  on 
the  three  sensor  probe  at  a  Reynolds  number  of  218. 


A 


Similar  tests  were  then  made  at  the  same  flow  speed  using 
smaller,  50.8  yU  m,  diameter  hot-film  sensors  (TSI-20).  The 
Reynolds  number  based  on  the  average  velocity  normal  to  the 
sensors  was  72.7.  The  variation  of  the  sensor  signal  with  yaw 
angle  was  found  to  be  more  regular  than  the  signal  from  the  152 
m  diameter  sensors*  but  sinusoidal  oscillations  of  the  sensor 
output  at  5.71  and  6.11  KHz  were  observed  at  small  yaw  angles. 
Even  smaller,  25.4yuV.rn,  diameter  hot-film  sensors  (TSI-10)  were 
then  installed  and  the  probe  was  again  tested  at  a  Reynolds 
number  based  on  the  average  normal  velocity  of  36.3.  The  sensor 
outputs  as  a  function  of  pitch  and  yaw  angle  were  much  more 
regular  than  the  152  yW.m  diameter  sensors,  but  6  KHz 
oscillations  of  the  sensor  output  and  jumps  in  the  mean  signal 


from  one  of  the  sensors  were  observed  at  large  pitch  angles. 


It  was  concluded  that  at  this  flow  speed,  the  Reynolds 
numbers  of  the  sensors  were  too  large  to  allow  steady, 
predictable  rates  of  heat  transfer.  The  aspect  ratio  of  the  hot- 
film  sensors  is  not  large  (of  the  order  of  15)  and  this,  as  well 
as  the  relatively  large  Reynolds  number,  probably  contributes  to 
the  irregularity  of  the  sensor  outputs  when  the  probe  is  pitched 
or  yawed.  It  should  also  be  noted  that  Ho13  has  reported  that 
electrical  cross-talk  occurs  between  the  sensors  on  a  split  film 
probe  with  diameter  of  152  jjjk.  The  present  finding  of  anamolous 
steady  and  unsteady  signals  resulting  from  a  change  in  the  angle 
between  the  flow  direction  and  the  sensor  axis  represents  an 
additional  source  of  aerodynamic  cross  talk  that  may  occur  with 
split  film  sensors  at  the  Reynolds  numbers  reported  above. 


III.  CALIBRATION  OP  A  MODIFIED  PROBE  PITTED  WITH  HOT-WIRE  SENSORS 


A.  Modifications  of  the  Probe 

The  TSI  1298  probe  was  slightly  modified  by  spreading  each 
of  the  three  pairs  of  hot-wire  prongs  further  apart  (in  a 
direction  outward  from  the  axis  of  the  probe),  as  shown  in  Fig. 
2,  to  reduce  wire/prong  interference  at  large  pitch  and  yaw 
angles.  The  approximate  diameter  of  a  circle  which  would 
surround  the  wires  at  this  larger  spacing  was  7  mm.  Thus, 
turbulent  velocity  fluctuations  with  a  spatial  scale  less  than  a 
centimeter  cannot  be  correctly  resolved  by  the  modified  probe. 
The  modified  probe  was  returned  to  the  manufacturer  and  fitted 
with  3.81  JlK.  m  diameter  tungsten  hot-wires  (TSI  T-1.5)  with  a 
platinum  coating,  a  sensitive  length  of  1270^tm  and  a  length  to 
diameter  ratio  333.3. 

B.  Determination  of  Probe  Geometry  and  Calibration 

Constants 

The  probe  was  tested  in  the  wind  tunnel  and  the  sensor 
signals  were  recorded  in  the  same  manner  as  discussed  for  the 
hot-film  probes  of  Section  II.  The  linearized  hot-wire  signals 
plotted  as  functions  of  pitch  and  yaw  angle  in  Figs.  5  and  6, 
were  much  more  regular  than  similar  traces  of  the  hot-film  sensor 
outputs,  see  Figs.  3  and  4.  The  signal  from  sensor  (2)  as  a 
function  of  varying  pitch  angle,  in  Fig.  5,  is  similar  to  the 
signal  from  sensor  (3)  as  a  function  of  yaw  angle,  in  Fig.  6, 
as  should  be  observed.  Further  examination  of  the  tape  recorded 
data  using  computer  generated  plots  of  the  linearized  output 


voltage  for  each  hot-wire  revealed  that  the  output  voltage  was 
proportional  (within  a  few  percent)  to  the  velocity  component 
normal  to  the  hot-wires  (the  classic  cosine  law).  The  output 
from  each  sensor  was  assumed  to  obey  the  cosine  law  so  that  each 
hot-wire  measures  the  component  of  the  velocity  normal  to  the 
wire.  For  three  sensor  probes  fitted  with  lower  aspect  ratio 
sensors,  which  do  not  obey  the  cosine  law,  Paulsen7  and  also 
Chang  and  Adrian®  used  the  Jorgensen®  cooling  law  to  relate  the 
sensor  output  to  the  flow  velocity  components. 

Using  the  geometrical  arrangement  of  the  hot-wires  shown  in 
Fig.  2,  The  equations  for  the  velocity  components,  U^,  normal 
to  the  ith  hot-wire  in  terms  of  the  velocity  components  relative 
to  the  probe  are, 

Ux  =  (  u2  +  v2  )1/2  (4) 

U2  =  C(  u  sin(p  -  v  cos  (p  )2  +  w2  31/2  (5) 

U3  =  [(  u  sin©  -  w  cos©  )2  +  v2  31/2  (6) 

/ 

where,  (p  ,  is  the  angle  between  sensor  (2)  and  the  probe 
axis,  ©  ,  is  the  angle  between  sensor  (3)  and  the  probe  axis 
and  u,  v  and  w  are  the  velocity  components  parallel  and 
normal  to  the  probe  axis  as  defined  in  Fig.  1. 


The  linearized  output  voltage,  ,  from  the  ith  hot-wire, 
is  proportional  to  the  velocity  normal  to  the  hot-wire,  i.e.: 

vi  =  Kiui»  (?) 

where,  Kj.#  are  constants. 

The  constants,  K^#  determined  from  the  digitized  records 
of  the  linearized  outputs  for  the  probe,  when  pitched  at  zero 
yaw  or  yawed  at  zero  pitch,  were  computed  by  noting  the  maximum 
output  voltage  for  each  wire  occuring  when  the  wire  was  normal  to 
the  flow.  Close  examination  of  the  digitized  data  revealed  that 
the  flow  in  the  tunnel  was  not  exactly  uniform  and  the  probe 
geometry  was  not  exact.  Hot-wires  (2)  and  (3)  did  not  appear 
to  be  mounted  accurately  at  45  degrees  to  the  probe  axis  or  else 
the  potentiometers  used  to  determine  the  angles,  ,  and,  ©y  , 
were  not  correctly  oriented  on  the  servo  shafts. 

An  iterative  procedure  was  developed  to  accurately  determine 
the  effective  probe  geometry  and  alignment  of  the  probe  axis  with 
the  flow.  Eqs.  (1-3)  were  used  to  calculate  the  velocity 
components,  u,  v,  and  w  as  functions  of  pitch  and  yaw 
angles,  Qp,  and  Qj  .  Then  Eqs.  (4-6)  were  used  to  compute  the 
velocities  normal  to  each  wire.  The  calculated  values  of 
velocity  normal  to  each  wire  must  agree  with  the  measured  values 
given  by  Eq.  (7)  using  the  previously  determined  constants,  K^. 
The  above  calculations  were  performed  using  various  assumed 
values  for  the  probe  geometry  (the  angles,  (p  ,  and  ©  )  and 

for  various  assumed  values  of  the  initial  alignment  of  the  probe 
axis  with  respect  to  the  tunnel  axis.  The  best  correspondence 
between  the  calculated  and  measured  normal  velocity  components. 


for  each  wire,  was  obtained  with  the  probe  geometry,  &  =  38 
degrees  and,  (p  -  45  degrees.  The  probe  yaw  axis  alignment 
ang 1 e ,  &y ,  was  found  to  be  correct,  but  the  pitch  axis 
alignment  angle,  0p ,  was  discovered  to  be  5  degrees  larger  than 
indicated  by  the  servo  pitch  axis  potentiometer.  The  above 
values  for  the  probe  geometry  and  alignment  of  the  probe  axis 
with  the  flow  were  used  in  the  remainder  of  the  paper. 

During  the  above  procedure,  the  response  of  sensor  (1), 
when  the  probe  was  pitched  at  zero  yaw  angle,  indicated  that  the 
mean  tunnel  speed  in  the  vertical  plane  of  the  test  section  was 
uniform  in  the  upper  part  of  the  test  section  (for  pitch 
angles,  0p  >  0°),  but  the  speed  was  higher  in  the  lower  half  (for 
a  pitch  angle,  ©p  =  -50°  the  speed  was  5%  greater  than  on  the 
centerline).  It  is  possible  that  the  flow  speed  varies  by  5%  at 
other  locations  in  the  test  section,  but  further  measurements  of 
the  flow  uniformity  were  not  performed.  For  this  reason,  the 
probe  calibration  is  not  as  accurate  as  it  would  have  been  in  a 
more  uniform  flow  and  the  measured  velocity  components  described 
below  may  be  in  error  by  5%. 


IV.  SIMPLIFICATION  OF  THE  THREE  SENSOR  PROBE  EQUATIONS 


Eqs.  (4-7)  must  be  solved  for  u,  v,  and,  w  when  the 
three  linearized  voltages,  V^,  are  measured.  From  Eqs.  (7)  and 
the  calibration  constants,  K^,  the  velocity  components  normal 
to  each  sensor,  U^»  are  known.  The  velocity  components,  u, 
v,  and,  w,  must  then  be  found  from  Eqs.  (4-6).  Numerical 


methods  were  used  to  solve  these  simultaneous  algebraic 
equations . 


The  equations  have  multiple  roots.  The  response  of  each 
hot-wire  depends  upon  the  magnitude  of  the  normal  velocity  and 
remains  the  same  when  the  velocity  is  reversed.  The  existence  of 
multiple  roots,  but  not  the  number  of  roots,  is  well  known  and 
has  been  discussed  by  many  investigators  including  Paulsen7  and 
Chang,  Adrian  and  Jones8.  Recently,  Willmarth14  presented  a 
geometric  interpretation  for  velocity  measurements  made  with 
three-sensor  probes  which  shows  that,  in  general,  there  are 
eight  possible  velocity  vectors  that  can  produce  a  single  set  of 
three  normal  velocity  signals. 

In  their  paper,  Chang,  Adrian  and  Jones8  used  the 
criterion  that  the  correct  velocity  vector  among  the  multiple 
results  is  the  one  with  the  least  angle  from  the  mean  flow 
direction.  They  found  this  to  be  valid  for  low  intensity 
turbulence,  but  for  highly  turbulent  flows  the  wrong  result  was 
obtained.  The  problem  has  not  been  resolved,  in  fact  Chang, 
Adrian,  and  Jones8  state  that,  "the  selection  of  the  right 
velocity  vector  among  the  multiple  results  becomes  a  major  task 
for  the  three  dimensional  velocity  measurement  in  higher 
intensity  turbulent  flow".  It  was  decided  that  the  best  solution 
possible  would  be  obtained  in  the  present  work  and  since  the 
correct  solution  is  known,  from  the  probe  orientation,  the 
limits  of  validity  of  the  solution  could  be  determined.  The 
problem  of  verifing  the  correctness  of  the  solution  in  a 
turbulent  flow  when  the  correct  solution  is  unknown  (as  in  the 
measurements  cited  by  Chang,  Adrian  and  Jones8)  is  left  for 


future  studies. 
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In  order  to  solve  Eqs.  (4-6)  for  each  set  of  the  three 
values  of  U^#  Eq.  (4)  was  first  solved  for  positive  values  of  u 


only. 


u  =  C  Ux2  -  v2]1/2. 


(which  eliminates  reversed  flow  solutions  along  the  probe  axis), 
and  substituted  into  Eq.  (6)  to  give  an  expression  for  w,  in 


terms  of  u,  and  v. 


w  *  u  tan 


9  -  C  U,2  -  v2]1/2sec  &  . 


The  sign  of  the  radical  in  Eq.  (9)  was  again  chosen  (in 
agreement  with  the  choice  of  sign  in  Eq.  8)  to  eliminate 
reversed  flow  along  the  probe  axis,  u  <  0,  when,  w,  and,  v, 
are  zero  or  small  compared  to,  U3.  Eqs.  (8)  and  (9)  were 
substituted  into  Eq.  (5)  to  obtain  an  expression,  depending  only 
on  v: 

F(v)=  C (U12-v2 ) */2sin <p  -  v  cos  <p  ]2 

+  C(U12-v2)1/2tan 0  -  (U32-  v2)1/2sec  &  ]2-  U22  =  0-  (10) 

When  Eq.  (10)  is  solved  for,  v,  the  components,  u,  and,  w, 
are  found  from  Eqs.  (8)  and  (9). 


V. 


ITERATIVE  SOLUTION  OF  THREE  SENSOR  PROBE  EQUATIONS 


Newton's  method  of  iterative  solution  was  used  to  find  the 
zeros  of  Eq.  (10).  For  an  initial  value  of  v,  say  vQ,  F(v0)  and 
dF(v©)/dv  are  used  to  find  a  better  approximation  of  v, 

vnew! 

vnew  -  vQ  -  F/(dF/dv).  (11) 

By  differentiating  Eq.  (10)  find, 

(dF/dv)/(2u)  ■  G[v/u,  (U12-v2)/(U32-v2)]  (12) 

where  the  equation  has  been  written  in  dimensionless  form  and, 

G(x,y)  ■  x[l  -  2sin 2  <p  -  (1+sin2  ©  )/cos2  ©  3 

+  x(  tan  0/cos  ©  )  (y  +  l/y)  +  (x2-  l)sin  (p  cos  (p  ,  (13) 

with:  x  *  v/u  and  y  »  [(U^2  “  v2)/(U^2  -  v2)]1/2. 

A)  Stability  of  the  iterative  solution 

Newton's  method  is  unstable  when  dF(v)/dv  *  0,  i.e.  G(x,y) 

*  0.  The  points  (w/u,  v/u)  for  which  G(x,y)  *  0  were  found 
numerically.  Since  G(x,y)  «  0  is  quadratic  in,  y,  for  any 
value  of,  x,  the 'two  roots  for  y  were  computed  numerically  for 
many  values  of,  x  ■  v/u.  Then  for  each  root,  y,  the 
corresponding  value  of  w/u,  was  computed  from  Eq.  (6), 
rewritten  (with  the  aid  of  Eq.  (4))  as, 

w/u  ■  tan©  -  y  sec©  (14) 

For  each  value  of  v/u,  there  are  two  values  of  w/u,  for  which 
dF/dv  »  0.  The  points  (w/u, v/u)  for  which  dF/dv  «  0  were 
computed  for  the  probe  geometry,  0  m  38°,  and  =  45°,  and 

plotted  in  Fig.  (7). 
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VI.  NUMERICAL  TEST  OF  THE  ITERATIVE  SOLUTION 


Two  FORTRAN  programs  were  written  to  test  the  performance  of 
the  probe  when  exposed  to  different  variations  of  the  flow 
velocity.  Numerical  data  for  the  normal  velocities  was  generated 
by  substituting  known  values  of  u,  v,  and  w  into  Eqs.  (4-6). 
These  data  are  guaranteed  to  be  correct,  unlike  experimental 
data  which  are  always  subject  to  error.  The  pecularities  of  the 
probe  geometry  and  the  iteration  procedure  can  best  be  examined 
using  numerical  data. 

Both  programs  first  calculated  the  velocity  components  u, 
v,  and,  w  using  Eqs.  (1-3),  then  the  velocity  normal  to  each 
hot-wire,  using  Eqs.  (4-6),  and  finally  performed  an  iteration 
procedure  to  calculate  u,  v,  and  w  from  the  normal 
velocities.  Each  iteration  was  carried  out  to  an  accuracy  in  v 
of  .01  ft/sec,  (approximately  0.02%  of  the  stream  velocity). 
The  starting  value  for  v  in  each  iteration  was  the  previous  value 
of  v  obtained  in  the  previous  iteration.  This  was  done  to 
increase  the  stability  of  the  Newton  iteration  procedure  and 
reduce  the  number  of  steps  in  each  iteration. 

The  first  test  program  alternately  incremented  and/or 
decremented  the  pitch  and  yaw  angles  in  one  degree  steps  creating 
a  "path"  in  the  pitch-yaw  plane  which  was  a  rectangular  "spiral" 
emanating  from  the  outward  from  the  origin.  The  spiral  continued 
outward  until  the  pitch  and  yaw  angles  reached  60  degrees.  The 
initial  starting  value  of  the  first  iteration  was  v  ■  0.  If  the 
iteration  failed  (negative  arguments  in  the  radicals  or  non¬ 
convergence)  the  iteration  was  restarted  at  the  next  point  on 
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the  spiral  with  the  initial  guess  for  the  value  of  v  being  the 
correct  value  of  v. 


The  u,  v,  and  w  solutions  found  by  iteration  were 
processed  to  eliminate  solutions  which  grossly  disagreed  with  the 
correct  values.  Fig.  8  shows  those  values  of  v/u  and  w/u 
corresponding  to  iterative  solutions  of  u,  v,  and  w  which  agree 
with  the  correct  values  to  within  5.0  ft/sec.  Interesting 
features  of  the  figure  include  the  rwo  triangular  shaped  holes  in 
the  first  and  third  quadrants  of  the  plots,  and  the  dark  band  of 
points  along  the  curved  upper  limit  of  the  fourth  quadrant. 

For  values  of  v/u  and  w/u  less  than  0.35  the  iteration 
procedure  was  reliable  and  the  resulting  solutions  were  quite 
accurate.  For  most  turbulence  measurements,  when  the  turbulent 
intensity  is  not  too  large,  these  results  are  more  than 
adequate . 

Closer  examination  of  Fig.  8  for  larger  values  of  v/u  and/or 
w/u  revealed  that  the  upper  slanted  edge  of  the  hole  in  the  first 
quadrant,  the  lower  slanted  edge  of  the  hole  in  the  third 
quadrant,  amd  the  ^urve  across  the  top  of  the  plot  in  the  fourth 
quadrant  are  all  segments  of  the  boundary,  dF(v)/dv  *  0  shown  in 
Fig.  7.  The  holes  in  the  first  and  third  quadrant  were  formed 
when  the  path  followed  during  variation  of  the  pitch  and  yaw 
angles  approached  the  boundary  from  the  side  opposite  the  origin 
and  the  solution  "reflected"  off  the  boundary  in  the  direction  of 
the  incorrect  solution  and  later  failed.  The  dark  band  across 
the  top  of  the  plot  in  the  fourth  quadrant  was  formed  when  the 
solution  path  approached  the  boundary  from  the  side  nearest  the 
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origin  and  reflected  back  toward  the  origin.  In  these  cases, 
the  solution  continued  quite  some  time  before  failing.  Since  the 
correct  solution  was  known  in  this  test  it  was  possible  to 
identify  the  errors,  but  it  would  be  very  difficult  to  determine 
the  correctness  of  the  solution  under  actual  conditions  when 
measuring  time  resolved  velocity  components  of  large  amplitude. 

A  more  demanding  test  of  the  iteration  procedure  was  made  in 
the  second  FORTRAN  test  program.  The  "path"  in  the  pitch-yaw 
plane  began  at  the  point  where  the  pitch  and  yaw  angles  are  -60 
degrees.  The  pitch  angle  was  incremented  in  one  degree  steps 
until  it  reached  60  degrees.  Then  the  yaw  angle  was  increased 
one  degree  and  the  pitch  angle  was  again  varied  from  -60  to  60 
degrees.  In  contrast  to  the  first  program,  if  the  iteration 
failed,  the  pitch  angle  was  increased  by  one  degree  and  the 
iteration  procedure  was  restarted  with  an  iniial  guess  for  v  of 


zero. 


Fig.  9  is  a  plot  of  the  values  of  v/u  and  w/u  generated  by 
the  second  test  program  which  were  within  5.0  ft/sec  of  the 
correct  values.  Again,  it  should  be  noted  that  the  iteration 
scheme  performs  quite  well  in  the  vicintiy  of  the  origin.  As 
expected,  most  of  the  good  points  lie  inside  the  boundary.  The 
use  of  the  initial  guess,  v  ■  0,  when  restarting  the  iteration 
after  failure  was  an  attempt  to  insure  that  the  solution  on  the 
side  of  the  boundary  nearest  the  origin  will  be  obtained.  This 
is  a  realistic  test  of  the  iteration  procedure  since  when  making 
actual  measurements  one  will  not  know  the  correct  solution  or 


where  it  is  located  relative  to  the  boundary. 


The  second  test  is  a  more  demanding  test  of  the  iteration 
procedure.  The  use  of  the  previous  value  of  v  as  the  initial 
guess  for  the  next  step  in  the  iteration  procedure  is  severly 
tested  when  the  pitch  angle  jumps  from  60  to  -60  degrees.  When 
restarting  the  procedure  with  the  initial  guess  of  v  *  0,  the 
iteration  procedure  was  observed  to  converge  to  solutions  lying 
on  the  side  of  the  boundary  dF(v)/dv  ■  0  adjacent  to  the  origin. 
A  comparison  of  Figs.  8  and  9  reveals  that  the  holes  in  the 
plots  do  not  coincide.  This  indicates  the  sensitivity  of  the 
iteration  to  the  starting  values  and  to  the  "path"  relative  to 
the  boundary.  A  test  was  made  in  which  the  iteration  procedure 
was  started  at  the  origin  and  the  yaw  angle  was  varied  from  zero 
to  negative  values  in  one  degree  steps.  The  correct  solution  was 
then  found  along  the  v/u  ■  0  axis  in  the  region  of  the  hole  in 
Fig.  9.  Apparently,  the  iteration  procedure  is  unstable  and/or 
sensitive  to  the  value  used  to  begin  the  procedure  in  this 
region. 


VII.  ITERATIVE  SOLUTION  FOR  MEASURED  DATA 


The  iteration  process  was  performed  using  hot-wire  signals 
measured  when  the  probe  was  manually  pitched  and  yawed  at  random 
for  a  constant  free  stream  velocity.  When  the  iteration  did  not 
converge  or  the  argument  of  one  of  the  radicals  in  Eq.  10  became 
negative  the  iteration  was  restarted  again  at  the  next  data  point 
with  an  initial  value  of  v  *  0.  The  probe  pitch  and  yaw  angles 
and  the  three  hot-wire  signals  were  recorded  on  analog  magnetic 
tape  and  digitized  off  line.  The  pitch  and  yaw  angles  for  this 
data  were  often  very  large,  greater  than  45°,  and  in  this  case 
the  iterative  solution  frequently  failed  to  give  a  solution  or 
gave  erroneous  results. 

Figs.  (10-12)  show  the  velocity  components  computed  by 
iteration  compared  to  the  actual  velocity  components  computed 
from  Eqs.  (1-3).  The  data  shown  in  Figs.  (10-12)  have  been 
screened,  using  the  computer,  to  eliminate  solutions  for  which 
any  velocity  component  differs  by  more  than  5  ft. /sec.  from  the 
correct  velocity  component  computed  from  the  probe  pitch  and  yaw 
angles,  see  Eqs.  (1-3).  The  data  of  Figs.  (10-12)  show  that  the 
probe  is  capable  of  measuring  velocity  components  produced  by 
flow  at  large  angles  to  the  probe  axis  with  reasonable  accuracy, 
but  it  is  essential  that  one  have  additional  data  to  reject 
incorrect  solutions. 

Another  test  of  the  accuracy  of  the  iterative  solution  was 
made  by  comparison  of  the  iteratively  computed  and  the  actual 
time  resolved  velocity  component  data  when  the  probe  was  randomly 
pitched  and  yawed  for  moderate  angles.  The  traces  shown  in  Fig. 


(13)  are  the  actual  and  the  iteratively  computed  velocity 
components  during  random  pitch  and  yaw  of  the  probe  at  a  wind 
tunnel  speed  of  approximately  65  ft. /sec.  The  velocity  traces 
have  been  plotted  so  that  the  origin  for  the  computed  velocity 
components  is  5  ft. /sec.  greater  than  the  origin  for  the  actual 
velocity  components.  The  iterative  solution  for  v/u  and,  w/u, 
is  also  shown  in  Fig.  14.  The  data  of  Fig.  14  do  not  reach  the 
boundary,  shown  in  Fig.  7,  where  dF(v)/dv  =  0  and  there  were  no 
failures  of  the  iteration  procedure  for  this  data.  However, 
some  of  the  data  are  incorrect. 

The  incorrect  data,  for  which  the  computed  and  true 
velocity  components  differ  by  more  than  5  ft. /sec.,  have  been 
replotted  and  flagged  (by  setting  the  iteratively  computed 
velocity  components  to  zero)  and  are  plotted  in  Fig.  15  with  the 
origin  of  the  correct  velocity  traces  5  ft. /sec.  above  the 
computed  velocity  traces.  The  streamwise  velocity  component,  u, 
was  fairly  accurate,  but  the  cross  stream  velocity  components 
were  in  error  at  large  amplitudes.  The  most  serious  error 
occurred  at  the  midpoint  of  the  data  where  the  true  and  computed, 
v  and  w,  traces  cross  over  each  other.  The  cause  of  the  errors 
are  discussed  in  the  next  section. 


VIII. 


DISCUSSION  AND  CONCLUSIONS 


When  the  probe  was  calibrated  for  this  investigation,  it 
was  found  that  the  stream  velocity  in  the  test  section  of  the 
tunnel  was  5%  greater  in  the  lower  portion  of  the  test  section 
than  in  the  upper  half.  This  was  not  deemed  significant  for  the 
measurements  shown  in  Fig.  13  and  15  because  the  pitch  and  yaw 
angles  of  the  probe  were  sufficiently  small  that  the  sensors 
remained  near  the  center  of  the  tunnel.  The  near  equality  of  the 
measured  and  correct  velocity  traces  for,  u,  shown  in  Fig.  13 
and  15  support  this  conclusion. 

However,  the  accuracy  is  not  as  good  for  velocity 
components  transverse  to  the  probe  axis.  When  either,  v/u,  or, 
w/u,  were  of  the  order  of  0.25,  or  greater,  a  large  number  of 
errors  were  found  for  positive  values  of,  v,  or,  w.  Some 
errors  were  also  found  for  negative  values  of,  v,  near  the  left 
side  of  Fig.  15  and  for  small  values  of,  v,  near  the  center  of 
the  trace  shown  in  Fig.  15.  There  appear  to  be  three  primary 
reasons  for  these  errors. 

1.  Instability  of  the  Method  of  Iterative  Solution 

The  instability  of  the  Newton  method  when  dF(v)/dv  =0  is  a 
serious  source  of  error.  Tests  of  the  method,  see  Sec.  IV, 
showed  that  the  iteration  placed  more  points  near  the  curves 
where  dF(v)/dv  =  0,  Fig.  8  and  9,  (hereafter  termed  the 
boundary)  than  elsewhere  on  the  plot.  However,  a  similar  plot, 
including  only  points  where  the  iteration  agreed  with  the  correct 
data  to  within  0.5  ft. /sec.,  showed  the  data  distribution  to  be 
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uniform  all  the  way  to  the  boundary.  The  extra  data  along  the 
boundaries  in  Fig.  8  and  9  were  found  to  correspond  to  values  of 
v/u  and  w/u  originally  lying  on  the  other  side  of  the  boundaries. 
This  shows  that  even  when  the  velocity  components  normal  to  the 
sensors  are  correct,  the  iteration  procedure  often  found  the 
incorrect  solution.  Observe  that  the  boundary  is  relatively 
close  to  the  origin  in  the  first  quadrant  of  the  w/u,  v/u 
plane,  passing  through  the  point,  v/u  =  w/u  =  0.38.  This  may 
be  why  most  of  the  errors  in  the  v  traces  and  all  of  the 
errors  in  the  w  traces  of  Fig.  15  occur  for  positive  values  of 
v  and  w. 

The  iteration  procedure  is  also  sensitive  to  the  initial 
starting  value  of  v.  The  test  in  Section  VI  could  not  find 
solutions  in  the  blank  part  of  the  third  quadrant  of  Fig.  8., 
but  the  correct  solution  was  found  when  the  iteration  was  started 
at  the  origin,  v  =  w  =  0.  The  erroneous  data  for  negative,  v, 
in  the  trace  at  the  left  hand  side  of  Fig.  15  may  have  been 
caused  by  an  instability  of  this  type. 

It  is  concluded  that  improved  accuracy  of  the  solution  is 

m 

possible  if  one  can  find  better  iteration  procedures  and/or 
methods  for  restarting  the  procedure  when  it  fails  to  converge. 
With  adequate  iteration  procedures  a  three  sensor  probe  can  be 


used  for  time  resolved  measurements  of  three  velocity  components 
in  a  turbulent  flow  if  the  turbulence  level  is  not  too  large. 
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2)  Incorrect  Measurement  of  the  Velocity  Normal  to  the 
Sensors 

For  certain  sets  of  velocity  components  normal  to  the 
sensors  which  correspond  to  points  (w/u,  v/u)  near  the  boundary 
there  may  be  no  real  solution  for  Eq.  (10),  if  one  or  more  of 
the  velocity  components  normal  to  the  sensor  is  erroneous. 
Errors  in  measurement  of  the  velocity  normal  to  the  sensors  can 
be  caused  by  errors  in  electrical  measurement,  aerodynamic  prong 
interference,  flow  nonuniformity  with  a  scale  less  than  the 
sensor  spacing,  poor  sensor  calibration,  dirt  on  a  sensor  and 
sensor  vibration,  see  Perry^. 

3)  Multiple  Solutions  of  Eqs.  (4-6) 

The  existence  of  multiple  solutions  for  a  given  set  of 
values  of  velocity  components  normal  to  the  sensors  can  result  in 
unavoidable  and  unrecognizeable  errors  in  the  solution  of  Eqs. 
(4-6).  Using  a  geometric  interpretation  for  measurements  made 
with  multiple  sensor  probes,  Willmarth14  found  that  there  are, 
in  general,  eight  velocity  vectors  which  can  cause  a  single  set 
of  three  velocity  components  normal  to  the  sensors.  Four  of 
these  solutions-  with  u  <  0,  have  been  eliminated  from 

consideration.  By  using  the  geometric  interpretation  one  can 
verify  that  there  are  always  two  possible  velocity  vectors  when 
the  transverse  velocity  components  are  small.  One  vector  is 
correct  and  the  other  is  an  incorrect  vector  at  a  large  angle 
with  respect  to  the  probe  axis.  One  can  then  choose  the  solution 

Q 

with  the  least  transverse  velocity.  Chang,  Adrian,  and  Jones 


recognized  this  for  their  probe  configuration  and  chose  the 
solution  with  a  flow  direction  making  the  least  angle  with  the 
probe  axis. 

When  the  transverse  velocity  is  larger,  there  are  either 
two  or  four  possible  solutions,  with  u>  0,  as  can  be  shown 
using  the  above  cited  geometric  interpretation  for  a  particular 
three  sensor  probe,  and  only  one  solution  is  correct.  In  this 
case  there  is  no  way  to  determine  the  correct  solution  of  Eqs. 
(4-6)  without  additional  information,  beyond  that  obtained  from 
measurements  made  with  three  sensors.  It  is  suggested  that  one 
might  use  additional  sensors  to  provide  more  information  so  that 
the  correct  solution  can  be  determined.  This  would  be  an 
interesting  problem  for  further  investigation,  but  would  make 
the  measurements  and  data  reduction  more  complex. 
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FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Line  drawing  of  servo  positioner  used  for  pitch  and  yaw 
calibration  of  the  three  sensor  probe  in  a  uniform  flow 
in  the  wind  tunnel.  The  housing  of  the  pitch  servo, 
(A),  is  fixed  and  the  housing  of  the  yaw  servo,  (B), 

is  mounted  on  the  rotating  shaft  of  the  pitch  servo. 

The  three  sensor  probe  is  in  turn  mounted  on  the 
rotating  shaft  of  the  yaw  servo  but  the  three  sensors 
at  the  head  of  the  probe  are  not  shown.  The  velocity 
components  with  respect  to  the  probe  are  (u,  v,  and 
w)  as  shown  in  the  sketch. 

Sketch  of  Modified  TSI  Type  1298  Three  Sensor  Probe 
Fitted  with  Hot-Wires  (TSI  model  T1.5).  The  velocity 
components  normal  to  each  sensor  are:  =  (u2+v  2)1/2, 

U2  =  [(u  sin  (p  —  v  cos  (p  )2+  w2]1/2,  and 

U3  =  C(u  sin©  —  w  cos©  )2+  v2]1/2 

Response  of  sensors  to  pitch  at  zero  yaw,  TSI  1298 
probe  fitted  with  hot-film  sensors,  (TSI-60).  Sensor 
diameter  =  0.152  mm,  length  =  2.032  mm.  Re^*  218. 

Response  of  sensors  to  yaw  at  zero  pitch,  TSI  1298 
probe  fitted  with  hot-film  sensors,  (TSI-60).  Sensor 
diameter  =  0.152  mm,  length  =  2.032  mm.  Re^*  218* 


Fig .  5 


Response  of  sensors  to  pitch  at  zero  yaw.  Modified  TSI 
1298  probe  fitted  with  hot-wire  sensors,  (TSI  T1.5). 


yt 


Sensor  diameter  =  0.0038  mm,  length  =  1.27  mm. 

Red=  5. 

Fig.  6  Response  of  sensors  to  yaw  at  zero  pitch.  Modified  TSI 
1298  probe  fitted  with  hot-wire  sensors,  (TSI  T1.5). 
Sensor  diameter  =  0.0038  mm,  length  =  1.27  mm. 

Red=  5. 

Fig.  7  Plot  of  values  of  (w/u,  v/u)  for  which  dF/dv  =*  0.  At 
these  points  the  iteration  scheme  is  unstable. 

Fig.  8  Results  of  test  of  the  iteration  procedure  when 
performed,  using  synthetic  data  for  the  velocity 
normal  to  each  hot-wire,  along  a  path  forming  a 
rectangular  "spiral"  emamating  outward  form  the  origin. 
Points  plotted  are  those  for  which  the  computed 
velocity  components  are  within  5  ft/sec  of  the  correct 
values. 
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Fig.  9  Results  of  test  of  the  iteration  procedure  when 
performed,  using  synthetic  data  for  the  velocity 
normal  to  each  hot-wire,  as  the  pitch  angle  was 
increased  from  negative  to  positive  values  at  constant 
yaw  angles.  Points  plotted  are  those  for  which  the 
iteratively  computed  velocity  components  were  within 
5  ft/sec  of  the  correct  values. 
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Fig.  10  Axial  velocity  component/  u,  computed  by  iteration 
compared  to  the  correct  value  obtained  from  Eqa.  (1-3). 
Pitch  and  yaw  angles  often  exceeded  45°,  thus  only 
iteratively  computed  values  within  5  ft/sec  of  correct 
values  are  shown.  Free  stream  velocity  »  65  ft/sec. 

Fig.  11  Same  as  Fig.l0f  for  the  transverse  component,  v. 

Fig.  12  Same  as  Fig.  10,  for  the  other  transverse  component,  w. 

Fig.  13  Time  resolved  traces  of  the  iteratively  computed 

velocity  components  compared  to  the  correct  velocity 
components.  Correct  velocity  component  traces  are 
plotted  displaced  5  ft/sec  above  computed  traces. 
Probe  pitch  and  yaw  angles  restricted  to  relatively 
small  angles.  Free  stream  velocity  approximately  65 
ft/sec. 

Fig.  14  Plot  of  iteratively  computed  values  of,  w/u  vs  v/u, 
for  data  shown  in  Fig.  13. 

Fig.  15  Same  data  as  shown  in  Fig.  13,  with  iteratively 


computed  velocity  component  differing  by  more  than  5 
ft/sec  from  the  correct  value  set  to  zero. 


Line  drawing  of  servo  positioner  used  for  pitch  and  yaw 
calibration  of  the  three  sensor  probe  In  a  uniform  flow 
in  the  wind  tunnel.  The  housing  of  the  pitch  servo, 
(A),  is  fixed  and  the  housing  of  the  yaw  servo,  (B), 
is  mounted  on  the  rotating  shaft  of  the  pitch  servo. 
The  three  sensor  probe  is  in  turn  mounted  on  the 
rotating  shaft  of  the  yaw  servo  but  the  three  sensors 
at  the  head  of  the  probe  are  not  shown.  The  velocity 
components  with  respect  to  the  probe  are  (u,  v,  and 
w)  as  shown  in  the  sketch. 


Fig.  2 


Sketch  of  Modified  TSI  Type  1298  Three  Sensor  Probe 
Fitted  with  Hot-Wires  (TSI  model  Tl.5).  The  velocity 
components  normal  to  each  sensor  are:  Ui  *  (u^+v2)1'2, 

U 2  *  £(u  sin  <p  v  cos  <P  )”+  w2]1/2,  and 
U3  *  C(u  sin©  w  cos  ©  )2+  v2]1'2 
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Fig.  4 


Response  of  sensors  to  yaw  at  zero  pitch,  TSI  1298 
probe  fitted  with  hot-film  sensors,  (TSI-60).  Sensor 
diameter  *»  0.152  mm,  length  ■  2.032  mm.  Red«  218. 
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Response  of  sensors  to  pitch  at  zero  yaw.  Modified  TSI 
1298  probe  fitted  with  hot-wire  sensors,  (TSI  T1.5). 
Sensor  diameter  -  0.0038  mm,  length  »  1.27  mm. 

Re<j«  5. 


Fig.  5 
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Fig.  6  Response  of  sensors  to  yaw  at  zero  pitch.  Modified  TSI 
1298  probe  fitted  with  hot-wire  sensors,  (TSI  T1.5). 
Sensor  diameter  =*  0.0038  mm,  length  =  1.27  mm. 

Red=  5. 
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Fig.  7  Plot  of  values  of  (w/u,  v/u)  for  which  dF/dv  =  0.  At 
these  points  the  iteration  scheme  is  unstable. 
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Fig.  10  Axial  velocity  component,  u,  computed  by  iteration 
compared  to  the  correct  value  obtained  from  Eqs.  (1-3). 
Pitch  and  yaw  angles  often  exceeded  45°,  thus  only 
iteratively  computed  values  within  5  ft/sec  of  correct 
values  are  shown.  Free  stream  velocity  =  65  ft/sec. 
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Same  data  as  shown  in  Pig.  13,  with  iterativel 
computed  velocity  component  differing  by  more  than 
ft/sec  from  the  correct  value  set  to  zero. 


